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making chain is fed by input parameters which are usually subject to uncertainties. The art of dealing 
with uncertainties has been developed in various directions and has recently become a focal point of 
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1. Introduction 


The uncertainty handling has been one of the main concerns of 
es the decision makers (including governors, engineers, managers, and 
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amraee@kntu.ac.ir (T. Amraee). energy sector decision makers are subject to a significant level of 
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e Line outage, 
e Generator outage 


Unemployment rates 
Inflation rate 

GDP 

Economic growth 


e Price levels 
e Governmental regulations 
e Market interdependencies 


Fig. 1. General classification of uncertain parameters in energy system studies. 


Fig. 2. Uncertainty modeling tools. 


There are various uncertainty handling methods developed for 
dealing with the aforementioned uncertain parameters as 
depicted in Fig. 2. The main difference between these methods is 
in line with the different techniques used for describing the 
uncertainty of input parameters. For example, fuzzy method uses 


data uncertainty [2]. The uncertain parameters in power system 
studies can be generally classified into two different categories 
including (see Fig. 1): 


© Technical parameters: These parameters are generally categor- 


ized in two main classes, namely topological parameters and 
operational parameters. The topological parameters are those 
related to network topologies like failure or forced outage of 
lines, generators or metering devices, etc. The operational 
parameters are tied with operating decisions like demand or 
generation values in power systems. 


membership functions for describing an uncertain parameter 
while the stochastic methods use probability density function. 
The similarity of them is that all of them try to quantify the effect 
of input parameters on model's outputs. These methods and the 
way the uncertainty is handled by them are described as follows: 


Economical parameters: The parameters which affect the © Probabilistic approach: One of the earliest works in stochastic 
economical indices fall in this category. Microeconomics inves- programming was done by Dantzig in 1955 [3]. It is assumed 
tigates the decisions of smaller business sectors like aggrega- that the input parameters of the model are random variables 
tors, domestic or industrial consumers while macroeconomics with a known probability density function (PDF). 

focuses on the entire power system industry. For example, © Possibilistic approach: The fuzzy arithmetic was introduced by 
uncertainty in fuel supply, costs of production, business taxes, Lotfi A. Zadeh in 1965 [4]. The input parameters of the model 
labor and raw materials are analyzed in microeconomics. On are described using the membership function (MF) of input 
the other hand, the issues like regulation or deregulation, parameters. 

environmental policies, economic growth, unemployment e Hybrid possibilistic-probabilistic approaches: Both random 
rates, gross domestic product (GDP) and interest rates are and possibilistic parameters are present in the model. 
analyzed in macroeconomics. All of these parameters are © Information gap decision theory: It was first proposed by Yakov 


subject to uncertainties and should be correctly addressed in 
economical studies. 


Ben-Haim [5] in 1980. In this method, no PDF or membership 
function is available for input parameters. It is based on the 


378 A. Soroudi, T. Amraee / Renewable and Sustainable Energy Reviews 28 (2013) 376-384 


difference between what is known and what is vital to be 
known by quantification of severe lack of information in the 
decision making process. 

© Robust optimization: It was first proposed by Soyster [6] in 
1973. The uncertainty sets are used for describing the uncer- 
tainty of input parameters. Using this technique, the obtained 
decisions remain optimal for the worst-case realization of the 
uncertain parameter within a given set. 

© Interval analysis: It was introduced by Ramon E. Moore in 1966 
[7]. It is assumed that the uncertain parameters are taking 
value from a known interval. It is somehow similar to the 
probabilistic modeling with a uniform PDF. This method finds 
the bounds of output variables. 


This paper is to provide a summary of recent techniques used 
for uncertainty modeling in power system applications. It offers 
a vision obtained from a relatively large number of previous 
works. This review serves as a road map to those interested in 
uncertainty modeling tools in power system studies to find the 
less explored research areas by standing on the shoulders of 
giants. 

The rest of this paper is set out as follows: Section 2 presents 
the Probabilistic approach, the possibilistic methodology is 
introduced in Section 3, the hybrid possibilistic-probabilistic 
approach is described in Section 4, the info-gap decision theory 
is explained in Section 5, the robust optimization technique is 
described in Section 6. Section 7 presents the interval analysis 
approach. Section 10 describes the promising lines of future 
researches. Finally, Section 11 summarizes the findings of 
this work. 


2. Probabilistic approach 


In the probabilistic approach, a multivariate function, namely y, 
y=f(Z) is available. Z is a vector of the form Z=[Z),...,Zm], in 
which z4 to Zm are random parameters with known PDFs while the 
PDF of y is tried to be identified. For a better explanation, the 
function f describes the system model (e.g. set of load flow 
equations), Z is a vector of input uncertain parameters to the 
system (e.g. power injections by renewable energy resources and 
electric loads) and y is the output variable (e.g. total active losses, 
total operating costs). Three probabilistic uncertainty modeling 
techniques are described as follows: 


2.1. Monte Carlo simulation (MCS) 


The Monte Carlo simulation is carried out in the following steps 
[8]. It is assumed that the z; are uncertain parameters. A sample, z?, 
is generated for each input parameter z;, using its PDF. The value of 
y® as the outcome variable, is calculated using y° =f(Z°) where 
Z? =[Z%,...,Z%,]. The procedure is repeated for a number of itera- 
tions, Nyc. Finally, the outcomes are analyzed using statistic 
criteria, histograms, confidence intervals, etc. There are some 
methods for reducing the computational burden of MCS like Latin 
Hypercube Sampling (LHS) [9], sample-splitting approach [10] and 
fission and roulette method [11]. 


2.2. Point estimate method 
The point estimate method (PEM) acts based on the concept of 


moments of uncertain input parameters. In a problem with n 
uncertain parameters, the major steps are as follows [12]: 


Step1: Set E(Y) = 0, E(Y°) = 0 and k=1. 


Step2: Determine the locations and probabilities of concentra- 
tions, €,; and Pki, respectively as follows: 


2 
1M i 1/M3(z 
Chi = 5 = HETI n4 3( - ») (1) 
k k 


Ek 3-i (2) 


2 
2ny/n+4 (e) 
2k 


where M3(z;) is the third moment of parameter Zp. 
Step3: Determine the concentration points z,;, as given below: 


1=1,2 (3) 


Py =(-1)! 


Zhi = Hz, t Eki X Ozi» 


where, 4, and oz, are mean and standard deviation of Zp, 


respectively. 
Step4: Calculate the f for both Zķi, as 
Z =i Zezen TH 1,2 (4) 


Step5: Calculate E(Y) and E(Y”) using: 


2 
E(Y) = E(Y)+ x Pkif (21,225 ---» Zki» ---»Zn) (5) 
2 
E(Y?) =E(Y?)-+ X Pif? Z1, Za» -<-s Zkis -< Zn) (6) 
i=1 


Step6: k=k+1 if k< n then go to Step 2; otherwise continue. 
Step7: Calculate the mean and the standard deviation as 


Hy =E(Y) (7) 


oy = \/ E(Y7)—E*(Y) (8) 


2.3. Scenario based decision making 


A scenario is defined as a probable realization of an uncertain 
set of parameters. A list of scenarios is generated using the PDF of 
each uncertain parameter, Z;. The expected value of output 
variable, y, is calculated as follows: 


y= È as x fZ) (9) 
seQy 


where }';<@,%s = 1 and zs is the probability of s scenario. 


If the number of scenarios are large then it is needed to obtain a 
small set of scenarios representing the original one. The purpose is 
to select a small set, Qs, with the cardinality of No,, from the 
original set, ©, [13]. A reasonable trade off must be respected 
between the loss of the information and decreasing the computa- 
tional burden [2]. The scenario reduction technique is carried out 
via the following steps [14,15]: 


Step1: Construct the probability distance matrix containing the 
distance between each pair of scenarios C(s, s) 
Step2: Select the fist scenario s4 as follows: 


s= arel min X messa) (10) 
SEOs EQ 
Qs ={S1},  Q)=Q)-Qs (11) 


Step3: Select the next scenario for Qs set, as follows: 


szar min X as min ass} (12) 


SE js e Qs) s" e QsU{S 
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Os = Qs U {Sn}, Qj = Q)-—Qs (13) 


Step4: If the cardinality of Qs is sufficient then go to step 2; else 
continue. 

Step5: Add the probability of each non-selected scenario to its 
closest scenario in the selected set, End. 


More detail can be found in [2]. 


3. Possibilistic approach 


Since the introduction of fuzzy set theory this technique has 
been used in many power system fields [16]. Suppose 
y=f(X1,...,Xn) is in hand and X vector contains the uncertain 
input parameters described using their associated membership 
functions. In this context, the function f describes the system 
model (e.g. self-scheduling problem for a genco in a liberalized 
electricity market), X is a vector of input uncertain parameters to 
the system (e.g. hourly electricity price) and y is the output 
variable (e.g. total profit of genco). 

Various membership functions can be used to formulate the 
degree of membership of a specific uncertain parameter depend- 
ing on the expert's opinion. Regardless of the membership 
function's shape the question is “how to determine the MF of 
y if MFs of X are known?”. The a—cut method can provide an 
answer to this question [17]. For a given fuzzy set A in U, the crisp 
set A% contains all individuals of U with membership value, A, 
greater than or equal to a, as calculated in 


A” = {x € U|u,(x) = @} 14 


A" = (AS A") (15) 


The a—cut of each uncertain parameter, x?, is determined using 
(14), then the a—cut of y, y%, is calculated as follows: 


y” =°,” (16) 
y” = (minf(X®), max f(X*) (17) 
X“ = (X°, X^) (18) 


In each a—cut, the upper bound of y, y”, and the lower bound of 
y%, y*, are maximized and minimized respectively. The final step is 
defuzzification. The process of translating a fuzzy number to a 
crisp one is called defuzzification [17]. Many defuzzification 
techniques are available such as maximum defuzzification techni- 
que, the centroid method [18], weighted average defuzzification 
technique. 


4. Hybrid possibilistic-probabilistic approach 


Sometimes, the decision maker is faced with a multivariate 
objective function, y=f(X,Z), where both possibilistic uncertain 
parameters (X) and probabilistic uncertain ones (Z) exist. 

For better clarification, the function f describes the system 
model (e.g. set of load flow equations), Z is a vector of input 
uncertain parameters to the system described by PDF (e.g. power 
injections by renewable energy resources and electric loads), X is a 
vector of input uncertain parameters to the system described by 
MF (e.g. electricity prices) and y is the output variable (e.g. total 
payments for procurement of active power losses). 

To deal with such cases some methods are developed which are 
described next. 


4.1. Possibilistic-Monte Carlo approach 


The following steps describe the mixed possibilistic-Monte 
Carlo approach [19]: 


Step1: For each z; e Z, generate a value using its PDF, zf. 
Step2: Calculate (y*)° and (y%)° as follows: 


(y%) = min f(Z°,X% (19) 
(y*)° = max f (Z°, X”) (20) 
X” = (X*,X") (21) 


These steps are repeated to obtain the statistical data of the 
parameters of the output's MF such as PDF or expected values. 


4.2. Possibilistic-scenario based approach 
The following steps describe this approach [20]: 
Step1: Generate the scenario set describing the behavior of Z, £2). 


Step2: Reduce the original scenario set to a small set, Q,. 
Step3: Calculate (y*) and Q”) as follows: 


y*=min ¥ as xJ ZX”) (22) 

= seQ, 

y*=max $ a; x f(Zs,X*) (23) 
SEQ; 

X” = (X%,X") (24) 


Step4: Deffuzzify the y. 
5. Information gap decision theory 


The Information Gap Decision Theory (IGDT) is a method to 
describe the uncertainties which cannot be described using PDF of 
MF due to the lack of sufficient information. It is used to make 
robust decisions against severe uncertainty of input parameters. 
Consider a typical optimization function as follows: 


y= min f(X, d) (25) 
H(X, d)=0 (26) 
G(X,d)>0 (27) 


where X is the vector of input parameters (which are subject to 
severe uncertainty) and d is the vector of decision variables. H and 
G are the equality and inequality constraints, respectively. f(X, d) 
describes the relations between the decision variables (d) and 
input uncertain parameters (X). 

In case the uncertain input parameters X are equal to their 
predicted values (X=X) then solving the (25)-(27) gives the 
predicted value of y=y. However, if the value of X is unknown 
then the IGDT method tries to find a solution for the problem 
which is robust against the error in predicting the value of X. In 
IGDT, the robustness is defined as the immunity of satisfaction of a 
predefined constraint [5]. The constraint satisfaction may be 
defined based on the application [21]. 

For a better clarification, assume that the function f describes 
the system model (e.g. set of constraints describing energy 
procurement from different resources), X is a vector of input 
uncertain parameters to the system which are subject to severe 
uncertainty (e.g. electricity price without any historic data) and y 
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is the output variable (e.g. total payments for energy procure- 
ment). d denotes the set of decision variables (e.g. amount of 
purchased energy from different energy resources like DG units, 
electricity pool market and bilateral contracts). The robustness in 
IGDT context is defined as follows: 

The total payments should be always less than a pre-specified 
threshold 7, no matter how the uncertain electricity price takes 
value far from what is predicted. The robust counterpart of the 
problem described in (25)-(27) is as follows: 


f(X, d) < če (28) 
fc=(14+0)xy (29) 
H(X,d) =0 (30) 
G(X, d) >0 (31) 


where ¢ is the degree that decision maker tolerates the deteriora- 
tion of objective function due to forecasting error of input para- 
meter X. The uncertainty of parameters in IGDT method is usually 
defined using the envelope bound model [22], as follows: 


X c U(a,X) (32) 


U(a,X) = Frise (33) 


where a is the uncertainty level of parameter X, X is the forecasted 
value of X and U(a,X) is the set of all values of X whose deviation 
from X will never be more than @X. The decision maker does not 
know the values of X and a. 

The robustness of a decision d based on the requirement /;, 
a(d,¢-), is defined as the maximum value of æ at which the 
decision maker is sure that the required constraints are always 
satisfied [5], as follows: 


åâ(d, lc) = max a 
S.t: Constraints (34) 


The decision making policy is defined as finding the decision 
variables, d, which maximizes the robustness, as 


max a(d, lc) (35) 
VX e U(a,X) (36) 
fX, d) < le (37) 
fc=(14+0)xy (38) 
H(X,d) =0 (39) 
G(X, d) >0 (40) 


6. Robust optimization 


The concept of robust optimization (RO) was first introduced by 
Soyster [6]. It is a new approach for solving optimization problems 
affected by uncertainty specially in the case of lack of full information 
on the nature of uncertainty [23]. It is described as follows: consider a 
function like z = f(X,y) which is linear in X and non-linear in y. The 
values of X are subject to uncertainty while the values of y are known. 
In robust optimization, it is assumed that no specified PDF is in hand 
for describing the uncertain parameter X. The uncertainty of X is 
modeled with an uncertainty set X e U(X), where U(X) is a set that 
parameter X can take value from it. The maximization of z = f (X, y) 


can be formulated via the following: 
max z =f (X.Y) (41) 


X € U(X) (42) 


Since the value of z is linear with respect to X, it can be reformulated 
as follows: 


max z (43) 
z<f(X,y) (44) 
F(X, y) =AW)*X +80) (45) 
X e U(X) = {X||X—X| < Å} (46) 


where X,X,X are the uncertain value, predicted value and max- 
imum possible deviation of variable X from X, respectively. The 
robust optimization seeks a solution which not only maximizes the 
objective function z but also insures the decision maker that if there 
exist some prediction error about the values of X, the z remains 
optimum with high probability [24]. To do this, a robust counterpart 
version of the problem is constructed and solved. The robust 
counterpart of (42) is defined as follows: 


max z (47) 
y 

Z<f(X.y) (48) 

F(X, y) =A aX +8(y)—Max Lai y)aX ew; (49) 

Lwi<T (50) 

0<w;<1 (51) 


Based on (47), two nested optimization problems are to be solved. 
Eqs. (49) and (50) are linear with respect to w; and has a dual form 
as follows: 


min p+ Dé) 
ŠP i 
P+ $i = aily)*x; (52) 
Inserting the (52) into (47) gives 
max Z (53) 
Y-i 
Z-Number 
e GY a iP 
lll” 
/ A 
/Fuzzy / Random\ 
/ Numbers 


Real Value Numbers 


Fig. 3. Uncertainty modeling trends: past, present and future. 
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Xj V <ô 


E <0 P+Qj 


Membership 
function 


PDF 


Universe of discourse 


Almost 


certainly Quit sure 


Membership 
function 


Ty 
Fig. 5. Concept of Z-number. 
z<fXy) (54) 
F(X, y) = AW)*X +80)-T p- Déi (55) 
B+E, = AW) ¥X; (56) 


There are some software developed for solving the robust optimiza- 
tion based problems [25]. As an illustrative example, consider that 
the function f describes the system model (e.g. set of constraints 
describing energy purchased by a smart home), X is a vector of input 
uncertain parameters to the system (e.g. electricity price which are 
always within a band), z is the output variable (e.g. total payments 
for energy procurement). Z” denotes the degree of conservativeness 
and y is the set of decision variables (e.g. amount of purchased 
energy in different hours). 


7. Interval analysis 


In this method, the range of values for each uncertain input 
parameter is defined and it can be represented by an interval. 
Suppose a multivariate function of the form f = (x1, ...,Xn) and 
lb; < xi < ub; where lb; and ub; are the lower and upper bounds of 
uncertain parameter x;, respectively. The goal is finding the lower 
and upper bounds of the objective function f. There are some 
softwares developed for solving the interval analysis based 


problems [26]: 


d 

1 292 
Prob= f A eT OH /20 
a o/27 
1 


b psi č dy 
_ X-a op ' f en")? /20? 4 | x d (xp? /20? 
oV2n\Ja b—a b c c-d 


(57) 


G(Prob) = pg, (Prob) (58) 


8. Exploring the new uncertainty handling methods 


The taxonomy of the uncertainty modeling methods in the 
past, present and future is as depicted in Fig. 3. In 2011, Zadeh 
introduced a new class of uncertain numbers called “Z-numbers” 
[27]. The Z-numbers are expressed as a pair in the form of 
Z=(A,B), in which, A,B are restrictions describing the behavior 
of Z. A is usually a fuzzy set while B describes the certainty degree. 
The certainty degree may be expressed as a PDF or a fuzzy set. In 
this context, Z = {x|x e A with certainty degree equal to B}. In clas- 
sic fuzzy numbers decision maker just has A and it is quit sure that 
Z belongs to A. However in Z-numbers, Z is described using the set 
A with a certainty (reliability) degree of information called B. 
Examples for Z-numbers are provided in Table 3. 

The normal PDF, as a function of jz, ø, is a reasonable choice for 
modeling the randomness of the load variable. In order to 
disambiguate this concept, a simple two-bus network is used as 
shown in Fig. 4. The series reactance of the transmission line 
connecting these two buses is assumed to be X. The voltage 
magnitudes of sending and receiving ends are represented by E 
and V, respectively. The angle by which the sending end voltage 
leads the receiving end voltage is considered to be 6. P and Q 
represent the active and reactive power load at the receiving end, 
respectively. The parameters of Z-numbers (pairs of MF) describ- 
ing load values are given in Table 4. 

For example, we are almost certain (set B2) that the demand 
value in a given bus (Z number) is low (set A;) as depicted in 
Fig. 5. The probability that the load value is low can be calculated 
as (58). 

In (58), G(Prob) indicates the degree to which Prob belongs to Aj. 
Now, the information of Z-number expressed as L = (A4, Bz) for load 
parameter is represented as a possibility distribution (G(Prob)) over 
the space of probability distributions (various values of p, o). 


9. Applications 


Context serves to demonstrate the applications of the afore- 
mentioned uncertainty modeling techniques. The applications 
are widely categorized into several fields, as given in Table 1. 
The summaries of uncertainty modeling attributes are provided in 
Table 2. 


© Distributed Generation (DG) impact assessment. 

è Plug-in hybrid electric vehicle (PHEV): (e.g. exploitation of plug 
in hybrid electric vehicles). 

e Assessment of available transfer capability (ATC). 

© Renewable energy (operation and planning) (e.g. hydro power 
generation management). 

© Load flow/optimal power flow calculations (e.g. probabilistic 
load flow, fuzzy load flow). 

è Reliability evaluation (e.g. reliability-oriented distribution net- 
work reconfiguration). 

© Distribution network operation and planning (e.g. phase bal- 
ancing, cost-benefit analysis of distribution automation). 
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Table 1 
Summaries of uncertainty modeling applications. 


Applications Probabilistic 
MC PEM 
DG units 28,29] 30 
PHEV 31,34] 35 
Available transfer capability (ATC) 37] 38 
Renewable energy (operation and planning) 40,41] 42 
Load flow/optimal power flow 50] 51] 
Reliability evaluation 56,57] 58 
Distribution operation and planning 29] 
Transmission/generation planning and operation/control 66] 
State estimation 78] 79 
Electricity market 84,85] 
* Unexplored research directions. 

Table 2 

Summaries of uncertainty modeling attributes. 
Method Input Output attributes Advantages 

representation 


Probabilistic PDF Statistics like expectation, variance, 
etc. 


Possibilistic MF MF 


Easy to implement 


Converting linguistic knowledge to 


Possibilistic Hybrid Interval RO IGDT 
Scenario 
[31,32] 33 [19,20] 2 
[35] 36] 
39 
[43-47] 48] 49] 
52-54] [55] 
[59,57] 60 [61,62] 
[63] 64] [65] 
[67,68,15,69] 70-73] [71] [74-76] 77) 
[80] 81,82] [81] [83] 
70 [86] 87-89] [90,88] 
Disadvantages 


Computationally expensive, needs a large amount of 
historic data, approximate result 
Complex implementation 


numerical values 


Useful for severe uncertainties 


Useful when just an interval is 


Dealing with both uncertainty types Computationally expensive 


Too conservative 


Controlled conservativeness 


The correlations among intervals are neglected this would 
make it too conservative 


© State estimation. 

e Electricity market (e.g. real time demand side management, 
bidding strategy, energy hub management and electricity 
procurement strategy). 

© Risk analysis (e.g. risk measures, risk hedging strategies). 


Hybrid MF & PDF Membership function with 
probabilistic parameters simultaneously 
IGDT Forecasted Decision variables satisfying the 
values requirements 
Robust optimization Intervals 
Useful when Difficult to 
justan use in non- 
interval linear 
is models 
available 
Interval Intervals Bounds of the outputs 
analysis available 
Table 3 
Examples for Z-numbers. 
Parameter A B 
Demand value High Very sure 
Wind speed Weibul PDF Normally 


Voltage magnitude Uniform distribution in [0.951.05] In most cases 


Table 4 
Describing the load values as Z-numbers. 


A B Load 
Low Not sure L = (Low, Notsure) = (A1, B1) 
Almost certain L = (Low, Almostcertain) = (A4, B2) 
Quit sure L = (Low, Quitsure) = (A4 , B3) 
Medium Not sure L= (Medium, Notsure) = (A2, By) 
Almost certain L= (Medium, Almostcertain) = (A2, B2) 
Quit sure L = (Medium, Quitsure) = (A2, B3) 
High Not sure L= (High, Notsure) = (A3, B1) 
Almost certain L = (High, Almostcertain) = (A3, B2) 
Quit sure L = (High, Quitsure) = (A3, B3) 


è Transmission/Generation planning, operation and control: (e.q. 
self-scheduling of GenCos, fault location scheme, dynamic 
economic dispatch, maintenance scheduling, determination of 
pilot points for zonal voltage control, small-signal stability). 


10. Promising lines of future researches 


The future trends in uncertainty modeling (to be investigated 
and further explored) are summarized as follows: 


10.1. Exploring new uncertain parameters 


With the increasingly revolutionary changes in power system's 
regulatory framework and developing technologies the uncer- 
tainty in input data of decision making procedures is increased. 
These uncertain environments include financial, societal/govern- 
mental (the ongoing government policy and the future potential 
incentive for the renewable energy), environmental (carbon emis- 
sion and global warming issue) and technical (communication and 
information architecture in smart grids, demand response, PHEV, 
energy hubs, smart building) uncertainties, risk preferences in the 
investment models, fuel prices and market regulations, renewable 
energy sources and competition among suppliers. 
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10.2. Enhancing the existing techniques 


© Reduce the computational burden specially when applied to 
large scale power systems and real-time applications. 

© Choosing the appropriate uncertainty handling technique. 

è Hybridizing the existing techniques to better describe the 
uncertain environment. 

© Using the heuristic methods to soften the computation procedures. 


11. Conclusion 


This paper proposed a standard classification of uncertainty 
handling methods along with the promising lines of future 
researches. The possibility of using Z-numbers for uncertainty 
modeling of load values was introduced for the first time. The 
assessed methodologies include probabilistic, possibilistic, hybrid 
methods, robust optimization, interval based analysis as well as Z- 
numbers. These models are compared and their strength and 
shortcomings are investigated. Based on the proposed comprehen- 
sive classification, it is deduced that each method is suitable for a 
specific type of uncertainty. The severity of uncertainty dictates 
choosing of the appropriate uncertainty modeling technique. 
Additionally, according to the carried out taxonomy of the meth- 
odologies, it was revealed that some research areas still remained 
untouched. 
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